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K. J. Arnold, M. P. Baden, M. D. Barrett 

Centre for Quantum Technologies and Department of Physics, 
National University of Singapore, 3 Science Drive 2, 117543 Singapore 
(Dated: March 4, 2013) 

We make a detailed experimental study of the threshold for self-organization of thermal 8T Rb atoms 
coupled to a high finesse cavity over a range of atom numbers and cavity detunings. We investigate the 
differences between probing with a traveling wave and a retroreflected lattice. These two scenarios lead 
to qualitatively different behavior in terms of the response of the system as a function of cavity detuning 
with respect to the probe. In both cases we confirm a iV -1 scaling of the threshold with atom number. 



Atoms coupled to the standing wave mode of a cavity 
will self-organize for sufficiently strong transverse pump- 
ing HI. In essence, light scattered into the cavity results in 
a potential that localizes the atoms to a configuration which 
favorably enhances collective scattering. Thus above a 
threshold pump intensity, an initially uniform distribution 
of atoms will undergo a phase transition, spontaneously re- 
organizing into a lattice configuration. Self-organization 
was first observed for thermal atoms in the experiments of 
ElO, and later with a Bose-Einstein condensate where it 
was mapped to the Dicke model 00. Self-organization 
is of particular interest as a platform for cooling as it can 
be applied to all polarizable particles, including molecules 
[3. In particular, theoretical studies have suggested 
cooling which has a rate independent of particle number 
N HI O, in contrast to other ensemble schemes where 
the cooling rate decreases linearly with N @[l0l. How- 
ever, numerical simulations suggested that the threshold 
may scale as N~ 2 instead of TV -1 if statistical fluctuations 
are required to trigger the self-organization JTll . For large 
ensembles, an threshold scaling places prohibitively 
severe constraints on the required probe power. Which 
threshold scaling applies will therefore greatly impact on 
the viability of self-organization as a cooling method lfT2l . 

In this Letter, we present a systematic experimental 
study of the self-organization threshold for 87 Rb atoms 
trapped in a high finesse optical cavity. We directly mea- 
sure the threshold behavior over a wide range of exper- 
imental parameters for two transverse probing configu- 
rations: a retroreflected lattice and a traveling wave, as 
shown schematically in Figure [T] In both configurations, 
the atoms are trapped intra-cavity by a 1560 nm far-off- 
resonance optical trap (FORT) which locates the atoms at 
every second anti-node of the 780 nm cavity mode. As 
discussed in Ref. [13], the threshold behavior depends on 
both the trapping configuration and the probing geometry 
used, resulting in a modification to the threshold equation 
in Ref. ifTTTl . However, the modified equations for both 
cases still maintain the TV -1 scaling within the mean field 
limit, which our experimental results clearly demonstrate. 

The threshold conditions relevant to our system are de- 
rived in Ref. [13] under the assumption that the system is 
in thermal equilibrium and that the transverse spatial ex- 
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FIG. 1 . Schematic representation of self-organization in the lat- 
tice (top) and traveling wave (bottom) geometry, (a, c) Below 
threshold, the atoms (black) are confined by the intra-cavity 1560 
nm FORT (yellow), (b, d) Above threshold, the atoms organize 
into a A-spaced lattice trapped by the probe and scattered fields 
(red), (b) For a lattice probe, the atoms can form one of two pos- 
sible A-spaced lattices (filled or open circles), (d) For a traveling 
wave probe, interference between probe and scattered fields re- 
sults in a A-spaced transverse lattice out of phase with the atoms 
by 0. 



tent of the atomic ensemble is large relative to the probe 
wavelength. For an atom-probe coupling fi, and detun- 
ing A = uj p — u a of the probe frequency uj p from 
the atomic resonance u a , threshold conditions are conve- 
niently expressed in terms of the dimensionless quantity 
H — —Titi 2 / (Ak B T), which is the probe trap depth rela- 
tive to the atoms' thermal energy, k B T. In our system, the 
threshold equations are given by fT3l 
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for the lattice geometry, and 
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for the traveling wave geometry. In these equations Ik(x) 
are modified Bessel functions of the first kind, ft is the cav- 
ity field decay rate, and A c = A c — NU^a. In the last 
expression A c = u p — uj c is the detuning of the probe fre- 
quency from the empty cavity resonance cj c , and NU^a 
is dispersive shift of the cavity resonance. The dispersive 
shift is due to N atoms each contributing a maximum sin- 
gle atom dispersive shift Uq — g 2 / A, where g is the atom- 
cavity coupling constant. The term a arises from averaging 
over the atomic spatial distribution. For a 1560 nm FORT 
with antinodes overlapped with every other antinode of the 
cavity mode, as in our system, a is given by fT3l 
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where 77 = V T o/(k B T), is the ratio of the trap depth, V T o, 
to the atoms' thermal energy. 

The experiments are carried out in a dual-coated high fi- 
nesse optical cavity. The cavity is 9.6 mm long and has a 
finesse F — 110, 000 near the wavelength of 780 nm and 
F — 160, 000 near 1560 nm. The high finesse at 1560 nm 
allows us to stabilize the length of the cavity as well as cre- 
ate a deep intra-cavity FORT. The intra-cavity FORT lat- 
tice has a waist of 70 /im and is actively stabilized to a 
trap depth of 230 fiK. The 1560 nm wavelength allows us 
to trap the atoms at exactly every second anti-node of the 
780 nm probe field, such that all trap sites are identically 
coupled to the cavity mode. The single atom cooperativity 
is C — # 2 /^7 ~ 6 as determined from the cavity QED pa- 
rameters (0,^,7) = 2tt x (1.1, 0.073, 3.0) MHz, where 
g is the single atom coupling constant for the | F = 2) 
to I F' — 3 ) cycling transition of the 87 Rb D2 line, ft is 
the cavity field decay rate, and 7 is the atomic dipole decay 
rate lfT4l . 

To load atoms into the cavity, we start from a magneto- 
optical trap (MOT) 15 mm above the cavity. We load up 
to 8 x 10 6 atoms into a single beam 1064 nm wavelength 
FORT similar to our previous experiments H21II6). This 
FORT beam is moved down 15 mm into the cavity over 
one second by a translation stage. Once in the cavity, the 
1064 nm FORT is adiabatically ramped off, transferring 
the atoms into the intra-cavity 1560 nm FORT. By varying 
the MOT atom number, we control the number of atoms 
delivered to the cavity FORT, up to a maximum of 7 x 10 5 . 

We probe the atoms with linearly polarized light which 
is aligned transverse to the cavity axis, as is a magnetic 
field determining the quantization axis. The experiments 
were performed at two probe detunings from the D2 tran- 
sition, — 110 GHz and —265 GHz. The atoms are opti- 
cally pumped into the | F — 1 ) ground state manifold 



resulting in a random distribution of | m ) states. How- 
ever with these detunings and polarization, the distribution 
across magnetic sub-levels is of no consequence. 

To verify the threshold equations, we need to measure 
the dispersive shift and the threshold intensity. The disper- 
sive shift is non-destructively measured by sweeping the 
frequency of a weak probe beam coupled to the cavity over 
the cavity resonance and observing the cavity transmission. 
Operating at a known laser-cavity detuning, A c , we are 
able to infer the dispersively shifted detuning, A c , from the 
directly measured dispersive shift, NU^a. Next the thresh- 
old is measured by linearly ramping up the intensity of a 
transverse probe beam over 10 ms while monitoring the 
cavity output. The cavity output is fiber coupled and split 
between a single photon counting module (SPCM) and het- 
erodyne detection setup. The SPCM is used for detecting 
weak signals, while the heterodyne detection is used to de- 
tect signals which would otherwise saturate the SPCM. The 
threshold for self-organization is clearly observed as a sud- 
den increase in cavity output, as illustrated in Figure |2ja). 

In order to properly compare the measured threshold 
power to the threshold equations three effects must be ac- 
counted for. First, the cavity has a significant birefringence 
separating the two linear polarization modes by 5ft. Lim- 
ited optical access constrains the probe beam polarization 
to be misaligned by 21° with respect to the cavity mode 
polarization. This reduces the scattering rate into the cav- 
ity by 13%. Second, in the lattice configuration, ramping 
up the probe beam results in a significant temperature in- 
crease. Since this temperature increase is not observed for 
the traveling wave probe, it is most likely due to adiabatic 
compression. As the threshold parameter fi depends on the 
temperature at threshold, we must calibrate this tempera- 
ture increase. This is done by measuring the temperature 
throughout an adiabatic ramp of the probe power with the 
cavity far detuned from the probe in order to avoid the on- 
set of self-organization. Third, the dispersive shift also de- 
pends on temperature via the parameter a. Thus, we com- 
pensate the measured dispersive shift using Eq[3]taking into 
account the initial temperature, 31 + 2 fiK, and the tem- 
perature at threshold as determined from the temperature 
calibration measurements. 

For the traveling wave geometry, Figure [2|c) shows the 
results of several threshold measurements at fixed probe 
detuning, A c . By coarsely controlling the mean disper- 
sive shift via the atom number, we sample a range of A c 
such as the data set shown. The threshold parameter is ob- 
tained by scaling the probe depth at which light is scat- 
tering into the cavity by the temperature at threshold. The 
black line shows the theoretical threshold calculated ab ini- 
tio from Eq[2] Data sets such as Figure [2jc) were taken at 
several values of A c ranging from — 65ft to— 5ft. In Fig- 
ure [2je) the black line shows the minimum threshold calcu- 
lated from Eq[2]at A c = and all measurements near the 
minimum for which | A c | < ft/2. As can be seen from the 
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FIG. 2. Threshold measurement results for the traveling wave 
and lattice probe geometries, in left and right columns, respec- 
tively, top: The self-organization threshold is determined from 
the sudden increase in cavity transmission as the probe intensity 
is ramped up over 10 ms, example traces are shown for a travel- 
ing wave (a) and lattice (b) probe, middle: One set of threshold 
measurements at fixed A c , for traveling wave (c) and lattice (d) 
probe. Black lines are an ab initio calculation of the threshold 
from Eq[T](Eq[2} for a lattice (traveling wave) probe, bottom: 
All measurements from the data sets at various A c near the mini- 
mum threshold, which for the traveling wave is where | A c | < k/2 
(e), and for the lattice is where | A c + k\ < k/2 (f). Black lines 
are the calculated minimum threshold for the respective cases. In 
(c-f) red circles are data taken at the atomic detuning —265 GHz, 
and blue circles at —110 GHz. 



figures, the threshold measurements are in good agreement 
with Eq.[5]and a N -1 scaling. 

For the lattice geometry, the probe beam is retro- 
reflected but we otherwise measure the threshold in the 
same way as for the traveling wave case. Figure^d) shows 
an example data set at fixed A c . The black lines are calcu- 
lated from Eq [T] where the minimum threshold occurs at 
A c = — k, in contrast to the traveling wave case. Fig- 
ure [2jf) shows all measurements near the minimum thresh- 
old, for which |A C + n\ < k/2. Again the results are in 
reasonable agreement with Eq[l]and a TV -1 scaling. 



However, in the lattice configuration there is an increas- 
ing discrepancy as the threshold parameter fi increases. We 
believe this to be a consequence of our trapping geome- 
try. Because the 1560 nm FORT restricts the axial mo- 
tion of the atoms, the probe lattice itself acts as a poten- 
tial barrier the atoms must overcome in order to organize. 
Clearly, in the limit that \i ^> 1 and the axial confine- 
ment is large, the atoms will not be able to organize at 
all. However, we observe a significant slow down in the 
transition to the organized phase near threshold even for 
/i approaching unity. We suspect the organization is ini- 
tially limited by the time required for the scattered poten- 
tial to reduce the barrier induced by the probe and facilitate 
runaway self-organization. This is supported by the obser- 
vation that organization again becomes rapid if we probe 
well above threshold. These effects result in the deviation 
of the measured threshold in Figure [2jf) as ji approaches 
unity. When we ramp the probe power, the slow onset of 
self-organization near to the threshold results in a delay be- 
tween crossing threshold and detecting a significant cavity 
output. Thus the threshold power is systemically overesti- 
mated. If the atoms were not confined in the axial direc- 
tion, we would expect rapid self-organization to occur for 
ji > 1. The absence of any systematic deviation in the 
traveling wave case, as seen in Figure |2|e), is consistent 
with this expectation. In this case there is no lattice poten- 
tial prior to organization so the atoms are not confined in 
the transverse direction. 

After the onset of self-organization, the cavity transmis- 
sion traces, for example Figure |5Ja-b), indicate complex 
dynamics for both probe geometries. In the case of the lat- 
tice probe configuration, our external potential geometry 
restricts us to ji < 1 for the reasons discussed above, and 
therefore we must operate in the regime of large dispersive 
shift (\NU a\ ^> k). In this regime a small change in 
NU a results in a significant change in A c . This results 
in non-linear dynamics above threshold due to the strong 
inter-dependence of the T, a, A c , and ji. We observe the 
organized phase persists for at most several milliseconds, 
as it does in Figure |2jb). We suspect an increasing tem- 
perature results in the system falling below threshold and 
reverting to the unorganized phase. This is supported by 
the data shown in Figure [3j Here the lattice was held at 
constant power above threshold and the temperature was 
measured via ballistic expansion after probing for 0.8 ms. 
The hysteresis in the temperature as the scattered poten- 
tial rises and falls indicates that the temperature rise is not 
simply due to adiabatic compression. It suggests the pres- 
ence of a strong heating mechanism possibly a result of 
non-adiabatic dynamics induced by the rapidly changing 
potential. 

In the case of the traveling wave, the dynamics above 
threshold are complicated by the fact that the induced or- 
ganizing potential is out of phase with atoms' positions by 
9 — tan _1 (— k/A c ), as discussed in Ref. lfT3l . When 
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FIG. 3. Indication of excess heating due to non-adiabatic dynam- 
ics after self-organization, (a) Measured temperature after self- 
organization by probing above threshold at constant power for 
0.8 ms plotted against the depth of the scattered field in the cav- 
ity as determined by the cavity output at the time of measurement, 
(b) Two example traces of the cavity output scaled the peak out- 
put corresponding to data points in (a). For red line (circles) the 
cavity output was above 50% of its peak at the measurement time, 
whereas for blue line (circles) the transmission was already below 
50% and falling. For reference, an intra-cavity photon number of 
4 x 10 6 corresponds to a 530 /jK deep optical lattice potential 
(A = -265 GHz). 



| A c | < k, this results in the organization switching off as 
the atoms are pushed away from the axial center of the cav- 
ity mode, decreasing the coupling to cavity and increasing 
the threshold condition. This effect can be seen from the 
cavity transmission trace Figure [2ja). Presumably as the 
atoms move back into the cavity mode, self-organization is 
re-initiated, resulting in the observed pulsing output. In 
principle the lattice can be stabilized with an additional 
force, as demonstrated for the CARL ifTTl . In our case, 
a displacement in the transverse Gaussian potential can 
counteract the phase and form a stable superradiant con- 
figuration. When | A c | ^> k, only small displacement is 
required since u ^ 0, and we indeed observe stable output 
for several milliseconds. 

In summary, we have shown that the self-organization 
threshold agrees well with a model based on simple mean 
field considerations given in Ref. jTSj . Earlier numerical 
work using multiparticle simulations indicated that scaling 
strongly depends on model details and assumptions ifTTIl . 
However we would argue that the mean field approach 
should yield accurate results for the threshold, provided 
there are sufficiently many atoms and a physical mecha- 
nism by which the atoms thermalize on the timescale of 
interest. In our experiments, the densities are such that col- 
lision times are ~ 100 /is with particle numbers of ^ 10 4 
per site of the 1560 nm lattice potential. In this case we 
can expect a thermodynamic description to be valid, and, 
indeed, the mean field results apply as evidenced by our 
measurements. For cases in which a thermodynamic de- 
scription may not apply, such as for fermions or low den- 
sities of molecules where the collision rate is negligible, 



a sub-iV -1 scaling of the effective threshold may still be 
applicable. 

Our experiments also indicated a strong heating mech- 
anism above threshold contrary to the cooling seen in the 
numerical simulations [11] but consistent with statements 
made in Ref. 0. In our case the heating is most likely due 
to the complex dynamics involved as the atoms organize. 
In future experiments, we will use a transverse 1560 nm 
lattice to confine the atoms. This will enable us to work 
more effectively in the regime of small dispersive shifts 
and to explore recent theoretical proposals for dissipation 
induced self-organization [8]. Additionally we will be able 
to establish a two dimensional A-lattice. In this config- 
uration, the atoms are expected to scatter super-radiantly 
without threshold into the cavity mode. In this context 
we can study potential cavity cooling mechanisms using 
weaker pumping and without the complex dynamics of the 
optomechanical forces driving the self-organized phase. 
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